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Abstract. The purpose of this article is to bring structure to (basic) hypergeometric biorthogonal systems, 
in particular to the q-Askey scheme of basic hypergeometric orthogonal polynomials. We aim to achieve 
this by looking at the limits as p — > of the elliptic hypergeometric biorthogonal functions from |13) , with 
parameters which depend in varying ways on p. As a result we get 38 systems of biorthogonal functions 
with for each system at least one explicit measure for the bilinear form. Amongst these we indeed recover 
the (j-Askey scheme. Each system consists of (basic hypergeometric) rational functions or polynomials. 



Elliptic hypergeometric functions are a generalization of hypergeometric series where the quotient of two 
subsequent terms is not a rational function of n (ordinary hypergeometric series), or a rational function 
of g" (basic hypergeometric series), but an elliptic function of n. They were introduced in the late 20th 
century, and since then many important identities for (basic and ordinary) hypergeometric functions have 
been generalized to the elliptic level. A nice overview of the study of elliptic hypergeometric functions is 
given in [14j . 

In [llj , [12) and |13| , Spiridonov and Zhedanov studied a family of biorthogonal elliptic hypergeometric 
functions. These functions can be seen as elliptic versions of the Askey- Wilson polynomials. The functions 
are defined as a very well poised elliptic hypergeometric sum 12^11 and the measure corresponds to the 
integrand of the elliptic beta integral. 

It should be noted that the elliptic hypergeometric functions do not represent polynomials, rather the 
biorthogonal functions are elliptic functions (of z) whose pole locations are restricted. Another significant 
difference is that on the elliptic level we consider biorthogonal functions, that is, we have two different families 
of functions, which arc orthogonal to each other. The two families of functions are related by interchanging 
two of their parameters. These two parameters do not occur for Askey- Wilson polynomials and determine the 
pole locations of the biorthogonal functions. For special values of the parameters a discrete biorthogonality 
exists, which can be obtained from the continuous measure by residue calculus. Just as the continuous 
orthogonality reduces to that of the Askey- Wilson polynomials, the discrete orthogonality reduces to that of 
the (7-Racah polynomials. 

The (7- Askey scheme 6 gives a description of all known univariate basic hypergeometric families of or- 
thogonal polynomials. It describes the families and the (limit) relations between them, and at the top we 
have the Askey- Wilson polynomials (i.e. all other families can be written as limits or special cases of the 
Askey- Wilson polynomials). In ^ the authors have shown that you can obtain a degeneration scheme of ba- 
sic hypergeometric beta integrals (both the ones with evaluations and those with transformations) by taking 
the limit as p Q from the elliptic hypergeometric beta integral and choosing how the parameters depend 
on p. While it is well known that a proper limit from the elliptic hypergeometric biorthogonal functions 
yields the Askey- Wilson polynomials, a systematic consideration of what limits can be obtained had not yet 
been carried out. Inspired by the results of our previous article we set out to discover what limits we could 
obtain, with the hope that the entire g- Askey scheme would appear as limits. 

We managed to obtain many different limits of (bi) orthogonal systems of functions, which includes the 
entire 5- Askey scheme. This is however just a subpart of our complete scheme. We also obtain numerous 
pairs of families of rational functions, which are biorthogonal to each other, families of rational functions 
biorthogonal to families of polynomials, and families of polynomials biorthogonal on the unit circle (e.g., the 
Pastro polynomials [8J). Several properties of the elliptic hypergeometric biorthogonal functions should be 
easily transferable (such as their behavior under difference operators), but we have not yet studied them all. 

One particularly pretty way of viewing all these systems and their relations is by considering the polytope 
P(0) from [2]. In that paper the authors tiled the polytope with three different tiles pf^ , P^f and . This 
tiling appears naturally in this paper by considering the Hesse polytope (P'^^ in [2J) restricted to a given 



hyperplane, and subsequently projected to another hyperplane. The tiles then correspond to the images of 
faces of the restricted Hesse polytope. 

As was shown in [2] each face of any of these tiles corresponds to some limit of the beta integral evaluation, 
which, as mentioned earlier is the squared norm of the constant function 1. Now we can extend this picture 
by associating to each simplicial face a system of biorthogonal functions (with orthogonality measure given 
by the appropriate beta integral limit). One system of biorthogonal functions is a limit of another one if the 
associated face contains the face of the other system in its boundary. When two faces differ by shifts along 
vectors in the root lattice of Eq, the associated biorthogonal functions are the same, though the two faces 
give two different measures. 

Let us now give a brief description of which limits we actually considered, and how we discerned in which 
cases the limits would form a biorthogonal system. We'll let Rn{z]tr;Ur',q',p) (where we have 4 parameters 
tr and 2 parameters Ur) denote the family of elliptic symmetric functions which is biorthogonal to the same 
functions with the two w-parameters interchanged. Now we consider limits as p 0, while the parameters z, 
tr and Ur of the biorthogonal function also change in this limit (but q remains fixed). In particular we look 
at Rn{zp'';trp"'';Urp'''';q;p). By a relatively simple argument, these limits always exist, for every choice of 

ar and ^r] that is, if we rescale the function by the proper power of p the limit as p — exists and is 
generically non-zero. However, many of these limits are independent of the variable z, which means that 
they cannot form a family of biorthogonal functions, thus we first filter those out. Our second criterion is 
to make sure that the limit of the bilinear form specialized to i?„ and Rm, i.e. the expression (Rn,Rm), 
still makes sense. Together these conditions determine our classification. It turns out that for each vector 
satisfying these two conditions we obtain a pair of families of functions which satisfy biorthogonality with 
respect to a measure we can write down explicitly. 

In [3] and |3] we show that the same classification also works multivariately, when we take limits of the 
multivariate i3C„-symmetric biorthogonal functions from [lOj and [9]. Apart from establishing orthogonality 
in those cases we also obtain generalizations of all other Macdonald conjectures there. The Macdonald 
polynomials will also fall in our scheme as special cases of multivariate Pastro polynomials. 

The article is organized as follows: First we have a notational section, followed by a general section on 
rings of power series in p, which is the space which contains the biorthogonal functions. Subsequently we 
discuss the definition of the elliptic biorthogonal functions in Section [3l Next we discuss the limits of the 
biorthogonal functions as p in Section |4l and in particular answer the question when these limits are 
z-dependent. In Section [3] we discuss systems of biorthogonal functions (two families of functions and their 
inner products) and consider when these have proper limits. Subsequently we briefly discuss the types of 
measures occurring in the limit. Section [7] compares our results with the g-Askey scheme, and the following 
section gives the relation with Pastro polynomials. Finally, in appendix |X] we explicitly tabulate all limiting 
systems. 



1. Notation of univariate ^-symbols 

We say a function f{x;z) is written multiplicatively in x if the presence of multiple parameters at the 
place of X indicates a product; and if ± symbols in those parameters also indicate a product over all possible 
combinations of + and — signs. For example 

n 

f{xi,X2, ■ . ■,xn;p) = Yi f{xt;p), 

r=l 

f{x^'^y^^]p) = f{xy;p)f{x/y;p)f{y/x;p)f{l/xy;p). 
Now we define the q-symbols and their elliptic analogues as in [5]. Let < |q|, \p\ < 1 and set 

oo n— 1 

{x]q) = Y[{l-xq''), {x;q)n^Y[{l~xq'^), {x;p,q) = JJ - xp^'ql 



r—O r,s>0 
n-1 



e(x\p) ^ {x,p/x\p), e(x\q\p)n = WO{xq'';p), T{x\p,q) = W_ j- 



r=0 i,j>0 



All these functions are written multiplicatively in x. Note that the terminating product (x; g)„ is also defined 
if |g| > 1. Likewise 9{x;q;p)n is defined for all q, though we must still insist on \p\ < 1. 

We call a function / : C* — )■ C symmetric if f{z) = f{l/z). A function is p-abelian if it satisfies 
f{z) = f{pz). 

We define the space A{uq]p^ q) as the space of all meromorphic symmetric p-abelian functions / such that 

+ l; ^ r^^^^) 



e{pqz /uo;q;p)rnff{z) ^ — +n /(^) 

1 [Uoq "^fz^'-) 



is holomorphic for sufficiently large mf (where my is allowed to depend on /). That is, / can only have 
poles at the points itog^'p*^ and u^^q^p^ for k £ Z and 1 < I < mf , and these poles must be simple. 

2. Function rings 

In this section we define some spaces of functions which contain the functions under consideration. By 
being slightly more restrictive in the space under consideration than usual we can prove some general results 
in advance. The variables x in this section will be specialized not just to the variable z of the biorthogonal 
functions, but also to the parameters tr and Ur- 

In principle we can consider the objects we work with as power series in a parameter p (with \p\ < 1). 
That is 

Definition 2.1. We let F = F{x) be the field of formal series J2teT '^ti^)P* ' where the set T of exponents 
is discrete and bounded from below, and the coefficients at{x) are (multivariate) rational functions in x 
independent of p, such that a^i^T ^ 0, together with the function which is constant 0. The valuation of 
a non-zero element f = '^^Otp^ ^ F is given by val(f) = minteTt. The leading coefficient is given by 
Icif) = ayaiij)- 

If the series f € F converges (for all small values of p) , then intuitively the valuation of / describes the 
asymptotic size of / as p — 0, while lc{f) describes limp_^oP^^"'^'^^/- This means that the goal in this 
article is to find leading coefficients (and to be certain we have the leading coefficient we need to find the 
valuations) of certain specific series (to be defined later). 

One of the properties we want to prove is the iterated limit property. Proposition [2?3l which is in essence 
a purely algebraic affair. In the statement of this proposition we must be able to substitute x by p'^x = 
{p'^^xi, . . . ,p'^"Xn)- To make sense of the new series f{p'^x) = X^teT 0't{p'^x)p* we must expand each at{p'^x) 
in terms of p and then add the resulting series together. This is not possible for every element in F{x); in 
particular the expansion of at{p'^x) may have an arbitrarily low valuation, which could give us infinite sums 
for the individual terms in the expansion of fip'^x), or violate the boundedness-from-below of the exponent 
set T. This problem can be avoided if we control the complexity of the rational functions at in an appropriate 
way. 

Definition 2.2. We define the degree of a rational function r £ C{xi, . . . , Xn), for any e €E M" as 

deg,(r(x)) = val{r{fx)), deg(r(a;)) = - M^^^^, 

where we use the supremum norm for ||e||. The subfield R{x) C F{x) is given by those series f = 
"^teT '^t{x)p'^ such that 

■ r deg(at) 

mf — — = G < oo. 

teT,t>vai{f) t — val(j) 

For monomials x" — Yii^^'^ ^'^e that degg(a;") = val{x'^p^''^) = e ■ u. Thus we find that deg(a:") = 
|itj|, so this is the usual definition of degree for polynomials, and we can now see that our definition of 
degree corresponds with (total) degree for polynomials. For univariate rational functions r{x) = x^^p{x)/q{x) 
(with a;, p and q being coprime polynomials and n G Z), the degree equals deg(r) = max(— n, n + deg(p) — 
deg(q)) (where we take the usual definition for degree of polynomials). Note that this degree can be negative, 
for example deg(a;/(l + x'^)) = —1, even though the degree of all (Laurent) polynomials is positive. It should 
also be observed that deg^{r{x) ■ q{x)) ~ deg^{r{x)) + deg^{q{x)) , so deg(r(a;) • q{x)) < deg(r(a;)) + deg((7(x)). 
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To see R{x) is indeed a field, there are two slightly tricky questions. First of all, addition of two elements 
in R(x) may reduce the valuation (if both leading terms have the same valuation), thereby bringing the 
denominator of deg{at)/{t — val{f)) closer to zero. However the discreteness of the exponent set T ensures 
that we still obtain a lower bound. 

The second tricky part is to show that we are allowed to divide in R{x). We can certainly divide in 
F{x), thus we only need to show the bound on the degrees of the rational functions. Indeed we have for 

/ \ 

^ p-val(f) p-val(f) 



and we have 



J]^ J^ptk-val(f) 



^ ^ Idj) 

ri>l (ti,...,t„)eT" fe=l ' 

tr>val{f) 



J 



LLi (^fe ~ '"o.lif)) k tk - val{f) T^iriteT.t>vai{f) t ~ a^(/) ' 



deg(ic(/)) 



In particular the constant C can increase by at most -. 7 rns- 

^ ■' miiitgr t>^„;(/) t~val{f) 

It can now easily be seen that for arbitrary / G R(x) and t ^ val(f) we have 

(1) valiat{p''x)p') ^ t + deg,{at) > < - ||e|| deg(at) > t ~ \\e\\{t ~ valif ))C. 

In particular if ||e||C < 1 the valuation becomes bounded from below, and in this case the series f{p'^x) is 
well defined. Note that we only know that f{p'^x) G R{x) for some, small enough e. In particular we are 
unable to consider f(p'^x) for arbitrary given e. 
Now the iterated limit property is almost trivial 

Proposition 2.3. Let f G R{x). Then for small enough e > and any u G R" we have 

lc{lc{f){p^x)) = /c(/(p™x)) 

and 

val{f) + e val{lc{f){p^x)) = «aZ(/(/"x)). 

Proof. This is clearly true for series f R consisting of a single term / = lc{f)p^°-^'--f\ If we add 
higher order terms atp*, they don't change the left hand side, and can only change the right hand side 
if val{at{p'^^x)p^) > val{lc{f){p'^^x)p'"°'^'^^^), as otherwise we have a unique lowest order term, determined by 
the leading coefficient. The explicit calculation ([T]) now shows that if e is small enough this is indeed the 
case. Indeed, we can choose e small enough such that t + \\4{t - val{f))C > val{lc{f){p''''x)p'"'^^f^) for all 
t eT with t> val{f). □ 

As a corollary we find that we can determine the valuation and leading coefficient of the sum of two terms 
in some cases where the valuations of the two summands are identical. Of course the case where the two 
valuations are unequal is trivial. 

Corollary 2.4. Let /, g G R{x) and define h = f + g. Suppose val{f) = val{g), and there exists a u ^ M" 
such that for all small enough e > we have val{f(p'^^x)) < val{g(p'^^x)) . Then val{h) = val{f) and 
lc{h)^lc{f) + lc{g). 

Proof. First note that as h{p'^"x) = + g{p'^"x) is the sum of two terms with different valuations, its 

valuation is equal to the minimum of those two valuations, i.e. val{h{p'^"x)) = val{f{p'^^x)). Now we can 
calculate using the above proposition 

val{h) + e • val{lc{h){p'^x)) = val{h{p'''^x)) = val{f {p"'' x)) = val{f) + e • val{lc{f){p'^ x)) . 

As this should hold for all small enough e > 0, we find that val{h) — val{f) (and val{lc{h){p'^ x)) ~ 
val{lc{f){p"x))). As a consequence we see that 

lc{h) = = + = + .gp-™'(f) = lc{f) + lc{g), 

where [p^] X^teT '^tP* ~ "^o denotes the coefficient of p^ in the expansion as power series in p. □ 



Of course we are not interested so much in power series, as we are in the (meromorphic) functions they 
define. Note that since we are allowing arbitrary real powers in the series, these are multi- valued functions 
of p. However, as in the case of left-bounded Laurent series, there is a radius R such that for < |p| < R, 
the series converges absolutely, and thus converges for any choice of branch. Thus the space we are really 
working in will be M(x), defined below. 

Definition 2.5. We say a formal power series f G F{x) converges uniformly on a compact set K , if 
p^'^'°-'-(i) f{^x) converges uniformly for x & K and p G B^{0) = {z \ \z\ < e} for some e > 0. M(x) is 
the intersection of R{x) with the field of formal power series f which converge uniformly to a holomorphic 
function for x in compacta outside the zero-set of some polynomial (which may depend on f). 

We also define A{uq]p, q) = A{uo]p, q) C\ M{uq) as the intersection of these functions with the symmetric 
p-abelian functions defined before. 

Note that the limit of a function / G F(x) is always its leading coefficient, and thus it is a rational 
function. This means that if / converges uniformly in a compact set K, than 1// converges uniformly on 
any compact set contained in K minus the zero-set of the numerator of its leading coefficient. In particular 
if / G M{x), then so is 1//, and therefore M{x) is indeed a field. 

Let us now consider some specific elements of M{x). First we define the constant function 

oo 

{p;pU = 11(1 -P') = ^(_l)«p"(3n-l)/2^ 
j=l nez 

which is an element of M where all rational functions (which are the coefficients of p*) are in fact ±1. Using 
Jacobi's triple product formula for the theta functions we find 

(p;p)oo^(x;p)-^(-x)>"("-i)/2^ 

nez 

so dividing these two elements shows that 6{x;p) G M{x) and l/6{x;p) G M{x). In fact we see that 
9{p°'x\p) G M{x) and l/9{p°'x;p) G M{x) for all a G R, by directly plugging in a; — t- p^x in the triple 
product identity. Note that for this implication we cannot plug in p^x after viewing 9{x\p) as an element of 
M(a;), as this could be ill-defined. We can easily determine the valuation and leading coefficient of 9{xp°';p) 
from this expression and obtain 

vamxp'';p)) = ^{a}{{a} - 1) - ^a{a - 1), 

where {a} = a — [aj denotes the fractional part of a, and 

lc{9{xp^;p)) 



(I - a G Z, 



We can now conclude that any sum of products and quotients of theta functions /, with arguments of the 
form p"a;" (for some monomial x'^) is indeed an element of M{x), and therefore has some valuation val{f) 
such that converges uniformly. 

3. BlORTHOGONAL FUNCTIONS 

In this section we recall the definition and basic properties of the biorthogonal functions from [13J. 
Definition 3.1. Let to, ti, t2, t^, uq, ui, q, and t be parameters such that tQtit2t3UoUi — pq. We define 

(2) Rn{z,to .ti,t2,t3,UQ,Ui,q,p) —l^^ftp.^.^^ gi-tpni qr^+Hp g +w. 

which is a very well poised elliptic hypergeometric series 12 Vn. 



The normalization for this definition is chosen so that the biorthogonal functions are highly invariant 
under shifts of the parameters. 
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Lemma 3.2. As functions of z we have Rn{z;ta : ti,t2,t^;uf),ui]q;p) G y4(?io;p,<7). Moreover the biorthogo- 
nal functions are elliptic in the tr and Ur, that is, they are invariant under multiplying these parameters with 
integer powers of p (as long as the balancing condition remains satisfied). Finally they satisfy the equations 

Rnizp^^'^-top^^^ : tip^'^^^,t2P^^^^,t3p^'^^'^;uop^^^,uip^^'^;q;p) = Rn{z]to : ti, is; "o, wii 

and 

~ 1 1 1 1 p p 1 ~ 
(3) Rn{z; — ; — , — ; = i?„(2;io : ti, ^2, is; uo, ui; 

to ti t2 ts Uo Ml q 

Proof. The biorthogonal functions are written as sums of functions in the space A("'(uo;p, q), so they are in 
this space themselves. The symmetries all follow from 6{px;p) = 9{^;p) = — i0(x;p). □ 

The biorthogonal functions satisfy a biorthogonality relation. There are two kinds of biorthogonality 
measures. For generic parameters we have a continuous biorthogonality, while for specific choices of the 
parameters the continuous measure reduces to a discrete one. The discrete version can be obtained from 
the continuous biorthogonality by residue calculus. The continuous version was discovered in [13], while the 
discrete version was already given in [TT] and [12] . While the biorthogonal functions themselves are perfectly 
well-defined for arbitrary q, the continuous measure below only works when |g| < 1 (otherwise the elliptic 
gamma functions in the integrand are not well-defined). The discrete measure is a completely algebraic 
affair, so it works for all values of q. 

Theorem 3.3. For any n,m E Z>o, and for generic values of the parameters such that io'^i^2^3'WoMi = pq 
we have 

{Rn{-;to '■ ti,t2,t3;uo,ui; q;p), Rm{-',to '■ ti,t2,t3;ui,uo; q;p))to,ti,t2,t3,uo,ui:q-p 

%,t2t3,^lt2,tlt3,^,^;g;p)„ 

'''"'Oi^;q;p)2n 9{^,tot„tot„tot2,j^^,j^;q;p)n'' 

where 

If \ f n \ ( ^ ^r^o^(^rZ^^P,g) dz 

{f.9hoMMMMM:,;v - 2Uo<r<s<,ntrt.;p,q) Jc^^"'^^"' r(z±2;p,,) 2mz' 

for parameters such that totit2t3tit5 = pq and functions f £ A(t4^;p,q) and g £ A{t5;p,q). Let mj be such 
that f(zi) T{tiZ^^) /T{tAq^™'i' zf'^) is holomorphic, and define rUg likewise for g. Let tr = tr for < r < 3 
and ti = t^q^^i and is = tc,q~™-!i . The contour is now taken to be symmetric (C = ) and contains all 
points of the form p^qHr (for i,j>0) (and hence excludes their reciprocals)^. 

If moreover toti = q^^ (implying no contour of the desired shape exists), and thus t2t'iUoUi = pq^'^^ we 
define the inner product as 

. - \- fu k.9{qtl;q;p)2k^itltotutot2,tot3,touo,i^;q;p)k ^ 

{hg)to,ti,t2,t3,Uo,Ui:q;p J [toQ )9[toq j" 



0{^:ht2,t,t3,'-^;q;p)N 
^ to ' ^^o^2 uot3 ' U(iui 

and have the same biorthogonality, unless q^ £ p^ for some < k < N (in which case one of the point 
masses becomes infinite). 

Note that using n — m = we find that the inner products are normalized such that (1,1) = 1. The 
corresponding equations are the famous elliptic beta integral evaluation, respectively the Frenkel-Turaev 
summation. It should be observed that one over the value of the inner products {Ri, Ri) gives the summand 
of a very well poised elliptic hypergeometric series (which, if it terminates can be summed by the Frenkel- 
Turaev summation) . This becomes convenient when studying determinacy of the measure after taking the 
limit to orthogonal polynomials. 



To be precise, C should be a chain representing the described homology class. 
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The definition gives us one expansion of i?„ as an elliptic hypergeometric series. It is known that the 
series in question, a 12VL1 has a large symmetry group, that is, that there are essentially 72 = |M^(i?6)|/|'S'6| 
ways of writing _R„ as such a series. Amongst these symmetries we find that the biorthogonal functions are 
symmetric up to a constant under permutations of the tr- The current normalization is such that 

(4) Rn{to;tQ : ti,t2,t3;uo,ui) = 1 

as setting z = to reduces the series to a single term. Together with the permutation symmetry in the tr this 
also gives us evaluations for Rn{tr) for r = 1, 2, 3. 

4. Limits of the biorthogonal functions 
The limits that we consider of the biorthogonal functions are the limits as p — >■ in 

for different values of ^, a^, 7r G R, where we assume tr, Ur and z are independent of p. That is, we explicitly 
describe how the parameters tr and Ur and the variable z depend on p. It should be noted that the balancing 
condition implies that + X^r 7i- = 1- Of course the limit will depend on the values of these exponents 

C, Qfr, 7r and we want to determine what the different possible limits are, and when those limits can give us a 
set of biorthogonal functions. The first necessary condition for obtaining biorthogonality in the limit is that 
the limiting functions must still be dependent on 2, so in this section we determine when this is the case. 

It should be observed that plugging in these values of the parameters in the definition of the biorthogonal 
functions, expresses i?„ as a sum of products of theta functions with arguments p^x for some x's and x's. 
By the general theory about power series we can thus immediately see that i?„ G M{z,tr,Ur,q), and thus 
that its limit exists (i.e. that it has some valuation and that 7?„p^"°''^^"^ converges uniformly to lc{Rn) as 
p outside the zero set of some polynomial) . 

What we do not know, however, is if we can plug in these values and then interchange sum and limit in 
the definition of the biorthogonal functions. Indeed there are many cases in which the valuation of i?„ is 
much higher than the valuation of the summands. Fortunately we have many different ways of writing the 
same function (using the symmetries of the 12V11), and in practical cases we can find a way of writing the 
function such that this is possible. 

For obtaining abstract information about the result of taking this limit, it is more convenient to express 
the series in terms of an elliptic beta integral. Indeed these series are equal to an from |2] specialized 
at a point where the integral reduces to a single series of residues. In that article the authors studied the 
limits of as p — >■ 0, so now we can apply the results from there. It turns out, we can obtain maximal 
symmetry in our equations by first using equation ([3]) of R„ which inverts q. This will impose the somewhat 
unconventional condition \q\ > 1, but since the results we obtain are about finite series, they are just as valid 
for arbitrary values of q. 

Therefore the expression in terms of i?^ we will work with is given by 



(5) Rniz; ta : ti, i2, ^3; uo, "i; 9, t; p) = 



pi( p^ p^q^ p^q^ ninh^v 1^ n^l 

q^"to q^'^tx g^"t2 q^ t3 ^ p£z 

' a ^ az^ a ^ OZ^ ' o"2: ' UnUi ' tnU-\ ' tnUri ' -f"^' ^ / 



q ^ qz^ q ^ qz ^ ' q^ z ^ UqUi ' taui ' taUQ 

■.ptl tot2 tpt'i . ( ^Og" P"0 *0"1 1 P£ 

q ^ q ^ q ^^^^ ' V Uq ' qtQ ' q^ ' UqZ^ Uq 



T-r3 p _p_ , , „_1 ptrZ U_ _Ptr_ -l\ (to. _PU0 , q^ Pq"Z IN 

In this expression we plug in z — > zp^, tr — > trP°^'' and Ur — >■ Urp'^^ ■ Subsequently the polytope P in 
the following proposition gives those vectors for which in [2 the limit of the above univariate biorthogonal 
function was given. In this proposition and the rest of this article we use the notation 04 := 70 and as := 71 
to indicate the symmetry between a^'s and 7r's. The symmetries of Lemma |3 . 2 1 change the vector (a; 7; C); 
For example ([3]) maps (a; 7; C,) 1—^ (— q;o, — ai, 1 — a2, 1 — as; —70, —71; C)- The proposition then shows that 
using these symmetries, any vector (a; 7; C) can be mapped into the polytope P. As a consequence we only 
have to consider the vectors (a; 7; C) ^ P and thus the limits for any univariate biorthogonal function can 
be found in [2]. 
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Proposition 4.1. The vector (i — ao, ^ — ai, . . . , ^ — as, ^ + C; ~^ ~ C) *^ ^'le polytope P^^'^ (from 'TJ 
if and only if (a; Q is in the polytope P defined by the hounding inequalities 



i 

Let T he the translation group with elements ta{x) = x + a, for a in the root lattice A of Eq, generated by 
the vectors 

^^i^i^iiiH^' (l'-l'0'0'0'0;0)' (0,1,-1,0,0,0;0), 

(0,0,1, -1,0,0;0), (0,0,0,1, -1,0;0), (0, 0, 0, 0, 1, -1; 0). 

Let G be the group generated by T and the flip 

flip{a; C) = (-ao, -"i, 1 - a2, 1 - as, -a^, -05; C)- 

Given any vector (a; C) G £ | X); '^'i — ^} ^ there exists an element g ^ G such that g(a; G P. 

Remark: Under the equation J^i Q^i = 1 the relation | > ai + aj + ak + \ ^ + C\ (for all < i < j < /c < 5) 
is equivalent to ^ + a,; + a^ + aj, > | i + C| (for all < i < j < fc < 5) and we will use them both. 

Proof. To obtain the bounding inequalities of the polytope P we just plug the variables (i— ao, ^— ai, . . . , i — 
as, 5 + C, — 5 — C) in the bounding inequalities of P^^^ and of the resulting equations we scratch those that 
follow from other equations in that list. The remaining equations are the ones we wrote down. 

To show that P contains a fundamental domain for G we first note that, using only the subgroup T, any 
vector V can be mapped to the polytope B defined by the bounding inequalities 

a, -aj<l, (l<i^j<6), a^ +aj > IC+ ^1 - ^, (1 < i < j < fc < 6), 

i 

Indeed by iterating the operation "add 1 to the smallest a^ and subtract 1 from the largest a^ if their 
difference is more than 1" we can ensure that a^ — aj < 1 (this operation must end at some point, as it 
decreases the non-negative integer J2i<jl\'^i ~ '^jW)- Now we note that (0; 1) is in the lattice A, so we can 
subtract or add an integer to ( to ensure |C + || < 1/2. If at this stage ai + aj + ak > \C + ^\ — ^ for 
all k, we are done. Otherwise the sum of the three smallest ai must be less than |C + ^| — ^- Suppose 
without loss of generality that ao > ai > a2 > as > a4 > as. Now consider the shift of a given by (/3; 9), 
where = a,. — ^ for < r < 2 and /J^ = a,, + ^ for 3 < r < 5 and 6* = C ± ^ , where we choose the sign 
to ensure |^ + ^| < ^ (which makes \( + ^\ + \6 + ^\ — We claim that now (/3; 9) G B. Indeed it is clear 
it still satisfies j3i — (3j < I (and + ^| < ^). Moreover we see that the smallest three /3's are exactly the 
largest a's minus ^, i.e. (32 and (3^. Moreover we have 

3 3 111 

l3o + I3i + 1^2 = ao + ai + a2 ~ - = I - as - a4 - a^ - - > -\C + - \ = \ - + 9\ - -. 

Now the polytope P is a subset of B, so we need to use the flip on the set B \ P and show that it 
lands in P. This is simplest if we use the ordered versions of P and B, that is we add the equations 
aa < ai < ■ ■ ■ < to the polytopes. We now write the complement as the union of two parts, the one 
part satisfying |^ + C| + ^ < ao + a4 + as , and the other part satisfying | ^ + CI + 5 ^ cto + Ck4 + Q^s- In 
the first case, in the complement B \ P, we can assume a4 + as > 1 (as if a4 + as < 1 and we were in the 
complement we would have ao < |^ + CI — 5 contradicting our new equation). Similarly, in the second case 
we can assume ao < |^ + C| — 5- 



Now we can look in the first case at the new parameters l3r — —ar (0 < r < 3) and l3r = 1 — ar {r — -i, 5) 
(and C = C), to land in the polytope with bounding equations 

/33</32</3i</3o</35</34<l + /33, 0<|i+C|<^, E^-^l' 

r 

/34 + /35<l, ^ >/3o + /3i+/32 + |C+^|, l>l3o + l34 + P5 + \C+l\- 

This is in the polytope P; the only equation which is not immediately clear is /Sa > || + C| — ^, which follows 
from 

/33 = IC+ ^1 - (/3o + A +/32 + IC+ ^1) + (1 -/34 -/^s) > IC+ ^1 - ^- 

In the second case we consider the new parameters /Sq = — ao and /3r — ^ ar for 1 < r < 5, and 9 
such that |C+^| + |0+^| = ^, and see that in that case too we land in P. 

We conclude that for any vectors in the complement of P in i? we can apply a flip and perhaps some 
more shifts to get a vector in P itself. □ 

So now we have to determine which limits are z-dependent. We first consider a slightly larger group of 
biorthogonal functions; those for which either the limit of is z-dependent, or the limit of the prefactor is 
z-dependent. 

Proposition 4.2. For a vector (a; Q G P, the limit of the in ([5]) depends on z for precisely one value 
o/ IC + 5I given a, to wit |C + ^1 = min(i, + 5, | — ctr — Q^s — a*) (where < r < s < t < 5). 

Proof. As shown in [2] the limit of E^ depends only on parameters orthogonal to the plane in which the 
face lies. The polytope in the parameters is a convex set, and for fixed values of a the possible values of C 
are given by some interval. Indeed if m = min(i, a^. + ^, ^ — Or — ag — at) (for all r < s < t) we find that 
C G [— m — i, TO — i]. The interval (— m — i, to — i) is now clearly contained in a single face of the polytope 
(as all inequalities with C are strict there), and this face contains the vector in the direction C. Hence the 
resulting limits are independent of z. Thus the only possible value of |C + ^| which gives a limit depending 
on z is IC + || = ™- n 

The z-dependent part of the prefactor can be given as 

{pz/uoA/uoz;p,q^^) 

(m"z/uo, q^/uoZ,puoz/q, uo/qz,puiz/q^, Ui/g"z;p, q~^) ]^r=a^P^rz/q, tr/qz;p, q-^) 

n-l 

= ]T L_ i 

11 (p^n fe2;/uo,g" V«oz;p) {puQz/q,Ualqz,puiz/q'\ui/q''z-p,q-^)Y{l^^{ptrz/q,tr/qz-p,q-^) 
We can now simply read off the following proposition: 

Proposition 4.3. Within the polytope P the prefactor above depends in the limit on z only if ar + 1/2 = 
IC + 1/2| for some < r < 5, or 04 -|- |C + i| > i. 

Of course, if both the prefactor and the Ei itself are z-dependent it is possible that their product is 
independent of z. Thus of the cases above we still need to check whether they are actually z-dependent or 
not. 

Proposition 4.4. Within the polytope P the limit of the biorthogonal function is z-dependent if and only if 
the vector {a; C) is in one of the following sets 

• The polytope P intersected with the half space defined by ai + + ^\ > ^; 



On the facet of P given by \C + ^\ = ^ + a4 + Ur + as for some r, s G {0, 1, 2, 3, 5}, r ^ 



2 1 2 

On the facet of P given by ar ^ \ = |C + ^| (for r G {0, l,2,3,5}j intersected with the halfspace 
04 + IC + ^1 1^ \j except in the interior of the resulting polytope; 
On the facet of P given by a4 + ^ — |C + 5I7 except in the interior of this facet; 



s; 



Proof. The limiting functions can only depend on z if either the Ei or the prefactor depends on z. This 
leads us to consider the subpolytope of P given by 04 + |C + 5I > ^, and the facets where |C + = 
IC + || = ar + (for some < r < 5), and \( + ^\ — ^ — a,- ~ as — at for some distinct < r, s,t < 5. 

To determine whether the functions in these polytopes are z-dependent, it suffices to show that one of 
their limits is z-dependent. We cut the polytope P according to the hyperplanes given by the equations 

+ ^1 = ^, IC + 5I = Or + 5, \C + ^ \ = ^ — ar — as — at and |C + 5I + 04 = 5- The iterated limit property 
shows that the limit associated to any vector (say in face F of P tiled using these hyperplanes) can be further 
degenerated to any face which has F as a subface. 

For example, take the subpolytope R which is given by the intersection of the half-space 04 + |C + 5I > 5 
with P. For the limit corresponding to any vector in R, we can always take a further limit to the limit in the 
interior of R. As the point (— |, |, |, §; |, |; —j) is in the interior, and as at this point only the prefactor 
depends on z, the function associated to the interior of R depends on z. Hence all functions associated to 
any vector in R are z-dependent. 

This also takes care of the case |C + = ^ and |C + ^| = | — — — 014 (with r ^ s <E {0, 1, 2, 3, 5}) 
as these facets of P are also facets of R. If we consider the facet of P given by |C + 5I = ^ + 04 + + Qfs 
(again with r ^ s ^ {0, 1, 2, 3, 5}), we note it is cut in two parts by the hyperplane 04 + |C + ^| = ^. On the 
part a4 + IC + 5I > \ we are in the subpolytope i?, whereas in the part 0:4 + |C + ^| < ^ gives a polytope 
(of codimension 1) in whose interior only the Ei depends on z. By again using the iterated limit argument, 
all limits associated to vectors on this facet of P depend on z. 

Let us now consider the case + | = |C + ^1 for some r ^ 4. This facet is again split into two parts by 
the hyperplane a4 + | + ^ | = ^, so by the same argument as above, we only need to consider the case for 
which 0:4 + + ^1 < 5. In the interior of the subpolytope of P given by these two extra equations, both Ex 
and the prefactor depend on z, and a direct calculation shows that their product in fact becomes constant 
(i.e. z-independent). As noted before, direct calculation of the limit of the Ex involves only plugging the 
parameters into the results of [2] (in this case using Proposition 4.3 in loc. cit.), and for the prefactor the 
limits are also immediate. So now we have to consider all facets of this polytope (which are codimension two 
in P), and direct verification shows that the limits associated to these facets all are dependent on z. Thus 
all limits in this polytope are z-dependent except the limits associated to the interior. 

The final case is when 04 + ^ = |(^-|- ^|. In this case we always have a4-|-|C-|-i| < i. A direct calculation 
(using [2| Prop. 4.3]) again shows that the limit of the biorthogonal functions is constant in the interior of 
this facet. Thus we calculate the limits of the facets of this facet and observe that those are z-dependent. 
Hence we conclude that the limit is z-dependent everywhere on this facet except on the interior. □ 

Finally obtaining the valuation of i?„ is simple, as we know the valuation of Ei to be zero (in the polytope 
P), and the valuation of the prefactor can be read off immediately. 

Proposition 4.5. For vectors (a; C) in the polytope P the valuation of _R„ is given by 
val{Rn{zp'^;top°'" : trp°''';urp'^'- j) = n(^{ao - a4)l{Qo-a4<o} - (-C - a4)l{-c-a4<o} 

- (1 + C - a4)l{l+C-Q4<0} - ("0 + a5)l{ao+a5<0} - ^{(^0 + ar)l{Qo-|-Q,<0} j 

It remains to figure out which limits have nice behavior in the inner product. 

5. Limits of biorthogonal systems 

We would like to consider the limits here of pairs of biorthogonal functions, together with their bilinear 
form. Of course the goal is to find new sets of biorthogonal functions, and a corresponding measure. 
For the measure we would like to have a limit equation of the form 

lc{{f {■;'trP°'"' ■,UrP^'-), g{-\trP°''' ■,UrP^'-))t^p<.^ ^Urp-<r) = {lc{f),lc{g))new,U,Ur., 

where lc{f) = lc{f{zp'^f ■,trP°'^;Urp'^^)) and lc{g) — lc{g{zp''i';trP°''';Urp'''')) for some C/ and (g. That is the 
leading coefficient of the bilinear form should only depend on the leading coefficients of the functions. Since 
the left hand side of this equation is symmetric in / and g, the right hand side must be as well, which implies 
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that wc must take Cf — Cgi that is, evaluate both functions around the same values of z. As the bilinear 
form is linear in scalars, it is obvious that the valuation of {f,g) depends on the valuations of / and g. But 
if those valuations are given the valuation of the inner product should no longer depend on / and g (at least 
for generic / and g). That is, we should have 

val{{f{-; trp"'';urp'''-), g{-; irP"''; WrP''''))t.p°'-,«.p^- ) 

= val{f{zp'^' ■,trp°''';urp'^'-)) + val{g{zp'^' lUp"'- ■,Urp'^'-)) + C(a,7), 
for some C{a, 7). But as we scaled our bilinear form to have (1, 1) = 1, we see that C{a, 7) = 0, thus we get 

(6) val{{f{-; Up"-'- ; Urp'^'' ),g{-; Up""- ; WrP'^'' ) ) t,p°r ,„,pTr ) 

= val{f{zp^i ; irP"" ; "rP'^")) + val{g{zp^i ; i^p"" ; UrP'^"))- 

It should be observed that we want an equation like this only for generic functions / and g. The case where 
/ and g are orthogonal would clearly violate such genericity (and in that case this equation would not make 
much sense). But expanding arbitrary functions in our biorthogonal families shows that taking / and g both 
equal to Rn (with interchanged parameters uq and ui) is generic enough. Thus we can consider the explicit 
biorthogonality norm relation 

{Rn{-'; to '■ ti,t2, ts] uq, Mi; q;p), i?,i(-; to : ti, t2, t^; ui, uq; q; p))to,ti,t2,t3,uo,ui;q;p 

^( ; 9; ^( ^ , toil ,tot,,tot2,J^^,J^;q;p)n'^ 

If we choose C (the location where we are evaluating the biorthogonal functions), then the valuation on 
the left hand side immediately follows from the valuations of the biorthogonal functions we obtained in the 
previous section. The valuation on the right hand side can easily be obtained. Thus this relation provides 
us with an equation between the valuations that must be satisfied in order for the system of biorthogonal 
functions to have a nice limit. The reason this equation might fail is if the valuation on the left hand side 
is not given by ([6|). With the usual rescaling of the bilinear form the limit of (i?,i,i?„) would then vanish. 
Thus in that case, either the limits of the biorthogonal functions are not linearly independent, or the limit 
of the bilinear form is not non-degenerate. Note that it can only be the case that the apparent valuation on 
the left hand side is less than its actual valuation, which equals the valuation on the right hand side, as only 
on the left hand side some extra cancellation can happen which we have not considered before. 

The valuation of the squared norm on the right hand side is a piecewise linear function of the a^, which is 
1-periodic in all (assuming the balancing condition remains satisfied), by ellipticness of the biorthogonal 
functions (and thus their norms). However it is rather complicated to write down explicitly for all ar- In 
the previous section we have seen that we could restrict ourselves to the case a G P (from Proposition 14. II) . 
so we will do the same here. In the polytope P the norm has valuation (remember 014 = 70 and =71) 

(7) - (70 +7l)l{7o+7l>0} - 2(70 +7l)l{7o+7l<0} + io'r + as)'i-{a^+a,<0} 
^ l<r<s<3 

3 1 . 

- ^(ar + ao)l{Q,-+Qo<0} + X! ^ "o)l{ao-7r>0} + ("O + 7r ) l{ao+7r >0} ) j • 

However, we are not so much interested in this valuation, as we are in this valuation minus the valuations of 
the biorthogonal functions which are the arguments of the inner product. The valuations of these biorthogonal 
functions are given in Proposition l4.5l Hence if we take the difference and divide by n, we find that we obtain 
an interesting limit if 

1 

(8) (70 +7l)l{7o+7l>0} + ("f- + "^)l{".+aa<0} "("^ + >)l{'rr+C>0} + (1 +C - 7r)l{l+C<7.} 

0<r<s<3 r=0 

equals zero. If this is positive the inner product will converge to zero (so we won't get anything interesting), 
and if this term is negative we must have made a mistake (as in that case a perfectly valid limit on the 
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left diverges on the right). Fortunately a (rather tedious) case analysis shows that latter case indeed never 
happens. 

It remains to determine when this term vanishes. One of the surprising results of the proposition below is 
that while both conditions (z-dependence and vanishing of this piecewise linear term) break the symmetry 
between the a^'s (0 < r < 3) and the 7r's, the final result is symmetric. 

Proposition 5.1. Consider the polytope P*-"-* (as in [2 J given by the bounding inequalities 
ar > — — , a,- — Q!s < 1, ar + as < I, a,. = 1. 

r 

Moreover define Pii,t (Pii.o = Pj^^ from 2 ) for < t < 5 as the polytope 

-^<at<0, at<ar<l + at, (r^t), < + < 1, (rT^sj^t), Qr = 1. 

r 

The set of vectors in P, for which each part of the associated limiting pair of the biorthogonal functions is z- 
dependent, and for which ([8]) vanishes consists of all vectors in P^*^) , outside the interiors of the subpolytopes 
Pii.t fO < t < with the value of ( given by \( + ^ \ = min(i, ar + ^,ar + as + at + ■^). 

Proof. Note that if 70 + 7i < then the sum ^ consists of no terms (and thus equals 0). Indeed we find for 
< r < s < 3 that ar + ag = 1 — at — a„ — (70+71) > (where t and u are such that {r,s,t,u} — {0,1,2,3}). 
And moreover 70 < —71 < | — |C+^| thus 70 + C < and 70 < C + 1 and similarly for 71. 

Moreover we see that at most three terms in the sum X]o<r<s<3 ' ' ' '^^^ non-zero, as an argument as 
before gives that either + ai > or 0:2 + 0^3 > (or both). This gives us 5 options up to symmetry: 

(1) + > for ah < r < s < 3; 

(2) ao + ai < and the sum of all other pairs is positive; 

(3) ao + ai < 0, ao + a2 < 0, and the sum of all other pairs is positive; 

(4) ao + ar < for 1 < r < 3 and the sum of all other pairs is positive; 

(5) ar + as < for < ?■ < s < 2 and the sum of all other pairs is positive. 

For the X]r=o ' ' ' P^^^ note that we cannot have both 70 + C > and 1 + C < 7i as 70 + 71 < 1. So 

there we have at most 2 non-zero terms and, up to symmetry (where we use ( ^ 1 — symmetry to choose 
-| < C < 0), 4 options: 

(A) 70,71 < -C,i + C; 

(B) 71 < -C < 70 < 1 + C; 

(C) -C < 70,71 <i + C; 

(D) 71 < -C < 1 + C < 70 

Within each of the 5 x 4 = 20 cases described by the above conditions the expression in ([5]) becomes linear, 
and thus reduces to a simple linear condition (though in case lA, i.e. the intersection of cases 1 and A, the 
condition is = 0). 

Now we recall from Proposition 14.41 that if the limit of the biorthogonal functions on both sides of the 
bilinear form is z-dependent and |C + i| ^ min(i, ar + i, ar -I- as -I- a* -I- i) (where we allow < r < 
s < t < 5) then we must have ^ — |C + ^| < 7o,7i, where we need both inequalities to ensure both 
elements of the pair are z-dependent. This puts us certainly in case C above. For each of the five different 
case for the signs of ar -|- as (0 < r < s < 3) the condition that the norm of the bilinear form has 
the right valuation gives us an equation for |C + i|. Inspection learns that in each case this equation is 
C -|- i| = min(i, ar -t- ^, ar + as + af -|- -i). To observe this we would like to remark that the inequalities in 
P imply IC + i| < min(i, ar + ^,ar + as + at + ^), so if, for example, in IC the equation ([5]) reduces to 
1^ -|- ^1 = ^, as a corollary we obtain that ^ — min(i, ar -I- ^ , ar -t- as -I- a* -I- ^). From now on we can thus 
assume that |C + ^| min(i, ar -I- ^, ar -f as + at + |) holds. 

It should now be observed that P*^°-' is the projection of P to the space of (ao, . . . , as). Thus now we have 
to show that for any vector in P(°\ except in the interior of one of the P//,t's we do satisfy all conditions, 
while in the interior of the P/z^t's one of the conditions fails. To show this we perform a case- by-case analysis. 
Our cases here depend on which of the following expressions (up to 6*4 x 5*2 symmetry) is minimal 

(i) h 

12 



(ii) ao + 5; 

(iii) 70 + 5; 

(iv) ao + ai + Q!2 + 5; 

(v) ao + ai + 70 + 5; 

(vi) ao + 70 + 71 + ^■ 

Note that Pn^t is exactly the polytope in which a* + ^ is minimal. Thus we want to check that for every 
vector in the facet of P given by term = |C + for term given by cases (i) and (iv)-(vi) above, the limits 
of the biorthogonal functions is z-dependent and the valuation of the bilinear form is correct (i.e. ([5]) 
vanishes). Moreover we want to check that outside the interior of the facets of P given by ao + ^ = |C + ^1 
or 7o + ^ = IC + the same holds, while in the interior either ([5]) does not vanish (case (ii)), or the limits 
of the biorthogonal functions are not z-dependent (case (iii)). 

In cases (i), (iv)-(vi) it turns out that we can identify which of the 20 (up to symmetry) linear pieces of 
(IH) we are (respectively IC, 5C, 2B and lA), and notice that the resulting linear equation is identically true. 
Moreover Proposition 14.41 shows that in these cases both biorthogonal functions have z-dependent limits, 
thus in these cases we are fine. 

In case (iii), corresponding to the polytope Pii,4, we are in the linear piece lA or IB of dS]), depending 
on the sign of 70 + 7i. In both cases we see that this linear function vanishes. However it is now obvious 
from Proposition 14.41 that in the interior of this polytope the first of the two biorthogonal functions has 
a z-independent limit (the second biorthogonal function has a z-dependent limit in this interior only if 
71 + IC + 5I > 5, or equivalently 70 + 71 > 0). Outside the interior the limits are all z-dependent. 

Finally case (ii), where ao + ^ is minimal and equals |C -I- ^| (i.e. we are in Pn^). First note that the 
functions are clearly z-dependent outside the interior of Qq. Let us now list all facets of P//.o up to the 
{Si X S'3) X S'2 symmetry on the vectors (ao, . . . , as). 

• ckQ = —1/2, we are in case 4A 

• — Q, we are in case IC 

• cko = ai, we are in case 2C 

• ao = 7o, we are in case IB 

• 1 -I- ao = as, we are in case 3 A 

• 1 -I- ao = 71 , we are in case 4B / 4D (in the intersection of these two) 

• ai + a2 = 0, we are in case 3C/5C (in the intersection of these two) 

• ai + 7o = 0, we are in case 2B 

• 7o + 7i = 0, we are in case lA 

• a2 + a3 = 1 , we are in case 2 A 

• as + 71 = 1, we are in case SB 

• 7o + 7i 1, we are in case 4C 

In all these cases we see that the valuation equation ([5]) linearizes and a direct calculation shows that it 
indeed vanishes. Thus outside the interior, all faces of P//,o correspond to a valid biorthogonal system. 

Finally we need to show that in the interior of P//,o the valuation equation ^ does not vanish. Thus we 
check all 20 polytopes on which (|8]) linearizes, add the equation ao + ^ — |C+ ^1 together with the vanishing 
of ([5]), and check whether these equations can be satisfied outside one of the facets of P//,o- To save us some 
work we note that in Pji^q we have ai + a2 > 0, so we can only be in case 5 if we are on a boundary of that 
case which is also covered in case 3. Moreover we note that 70 < ao + 1 = ^ + |C + 5I1 so we can also only 
be on the boundary of case D which is also covered in case B. Thus we only need to consider the 12 cases 
1234ABC. Apart from case lA the equation that ([5]) vanishes is exactly the defining equation of a facet of 
Pi I ft- In case lA on the other hand, we see that ([5]) vanishes as long as 70 + 71 < 0. As 70 + 71 > is a 
bounding equation for P//,o this implies that we must have 70 + 7i = 0, which again determines a facet of 
Pii.o- Thus we see that ([S]) vanishes nowhere in the interior of P/j.o- D 

We now know for which vectors the limit leads to two families of z-dependent functions with correct 
valuations so that we can hope that the limit from the elliptic hypergeometric level leads to a biorthogonality 
relation between the limiting functions. There are two more questions we need to answer. First of all we 
want to show that, indeed, for each of these vectors we do get a biorthogonal system. We do this by explicitly 
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defining the limiting measures and showing that we really can take a limit in the biorthogonality relation 
at the elliptic hypergeometric level. Secondly we want to know what different limits correspond to these 
vectors. As we saw in [2] taking limits along different vectors often give the same functions, so we want to 
tabulate what different limits we can obtain. 

The polytope P is just the intersection of the polytope P*^^-' with a given hyperplane (ay + as = 0), 
corresponding to the set of parameters in the Ei of (O, thus the limit of the Ei depends only on the 
face of P the vector (a; Q is contained in, and different faces of P lead to different limits of the Ei. The 
polytope P*^"-* is the projection of the polytope P by ignoring C. In this projection the inverse image of a 
face of P*^"^ correspond to a union of faces in P, and not always a single face of P. In particular, limits 
associated to vectors in the same face of P^") can still have different values. The point is that in the equation 

+ ^1 = min(i, ^ + ar, ^ + ar + as + at), C does not depend linearly on the vector a within each face 
of P^"^. However, if we restrict ourselves to the different facets of P given byC + ^ = ^, C+'^ — \ + ar 
(for < r < 5) and C + 5 = \ + ar + ag + at (for 0<r<s<i<5), and take the projections of these 
facets (called resp. P/, P//,t (the same as above) and P//7,(r,s,t)): the different faces in the projection will 
lead to different limits for the Ei . By inspection we can see that the limits of the prefactor in ([5]) also only 
depend on which face of one of these projections the vector is in. Therefore we see that the limit of the 
biorthogonal functions and their squared norm formula only depend on which face of P/, Pn^t and Piii^(r.s,t) 
contains the vector a (if we assume C, is given by the unique value for which we can have a biorthogonal 
system as limit). It could conceivably be that while the limits of the i?i's and of the prefactors is different 
in two different faces, the product of the limit of the Ei and the prefactor (i.e. the limit of the biorthogonal 
functions themselves) is the same on two different faces. By inspecting the different limits associated to the 
faces and considering basic properties of these limits (symmetries, locations of poles), we can exclude this 
possibility. Thus we obtain the following theorem. 

Theorem 5.2. Define the polytopes Pj, Pii,t (0 <t < 5) and Piii,(r,s,t) (0 < r < s < t < 5) as in ^ , that 
is 

• Pj is given by the bounding equations 

5 

ar > 0, ar = 1. 

r=0 

• Pii,t given by the bounding equations 

1 ^ 
--^<at<0, at<ar<l + at (r^t), < + < 1 (r, s/i), ^a^^l- 

• Piii^(r,s,t) is given by the bounding equations 

5 

aa + Q:6<0 (a, 6 e {r, s,t}), -{ar + as + at) < aa < I + {ar + as + at) {a^{r,s,t}), Qg = 1. 

a=0 

Then these polytopes tile P^*^-* . The leading coefficients of the triple 

{Ri{z;tr;uo,ui; z) , Ri{z;tr;ui,uo), {Ri{-;tr;uo,ui), Ri{-;tr;ui,uo))t^^u,.), 
when we replace tr — > trp"'' , Ur p°'T+i ^.f^^ ^ ^ pC y^lien C is determined as a function of a by + ^ = 
min(i, i + Q!r, 5 + + + at), for a vector a G P^^-* is determined by the face of the tiling of P*^°-' given 
above which contains a. 

If two faces are related by a shift from the translation group T from Proposition the corresponding 
triples have identical leading coefficients. If two faces are related by the action of x S2 on P^*'-' the leading 
coefficients are equal up to the same permutation of the parameters tr and Ur- If the faces are related by the 
flip of Proposition (followed by a translation in T ), the leading coefficients of the biorthogonal functions 
are related by setting tr — > l/tr, Ur — ^ 1/ur and q 1/q, in this case the squared norms are identical. 

We are left with the problem that we want to find (at least one) explicit bilinear form turning each of 
these cases into a biorthogonal system. It turns out that in [2] we determined limits of the elliptic beta 
integral evaluation (i.e. the equation (1, 1} = 1) for parameters specialized like this as p — )> 0. The kind of 
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limit depended precisely on the face the vector a was in in the tiling of P^") as in the above theorem. These 
limits all essentially involved evaluating the integrand at zp'' for z in some set which is independent of p 
and C such that IC + ^| = min(i, + ^,ar + as + at + that is, exactly where we want to evaluate our 
biorthogonal functions. It turns out that the exact same limits are still valid if we plug in the biorthogonal 
functions. Thus in all the cases mentioned in the theorem above we obtain at least an algebraic measure 
giving biorthogonality. In the next section we will describe these limiting measures briefly. 

In the cases of biorthogonal functions associated to a vector not contained in any translate of P^'^\ that 
is, the flips of 0022pp, 04as (i.e. Continuous (/"^-Hermite), and OOSlas (i.e. Stieltjes-Wiegert), we can 
obtain measures by taking limits outside this polytope. These measures are related to bilateral series. See 
[U Section 5] for details. 

It should be noted that the measures obtained are not necessarily positive. Of course you would only 
hope this were the case if you were considering orthogonal polynomials (i.e. cases where the limits of both 
families of biorthogonal functions are identical) and the squared norms were positive. However even in those 
cases the measures we find are not always positive. We believe you can find positive measures in those cases 
by taking limits from the elliptic level by looking at limits outside the polytope as in [31 Section 9]. 

Finally we should make a few remarks about explicitly obtaining the limits of the biorthogonal functions. 
Our description so far shows that we can obtain a limit by first rewriting the biorthogonal functions as beta 
integrals. The resulting limit will often again be a singular integral, which reduces to a finite sum of residues. 
This is far from our desired method of looking at the defining sum ([2]) for the biorthogonal functions and 
replacing sum and limit. Unfortunately the latter method does not always work (essentially because the 
valuation of the sum is more than the valuation of the summands). However due to the VF(i?7)-symmetry 
of the elliptic hypergeometric series there are many different ways of writing the biorthogonal functions as 
a series. In [3] we show that 8 of these representations generalize to the multivariate level, and at least 
one of those admits interchanging sum and limit. In particular, if we want to do explicit calculations we 
do not have to make the detour through singular beta integrals in order to get expressions for the limiting 
biorthogonal functions. 

Finding all different biorthogonal systems of rational hypergeometric functions which are limits of the 
elliptic hypergeometric biorthogonal functions has now been reduced to a combinatorial exercise of writing 
down the different faces of the given tiling of P^^^ modulo the symmetry group 5*4 x ^2. Each such face 
corresponds to a pair of limiting biorthogonal functions and a corresponding limit of the measure. If two 
faces are related by a translation from T in Proposition 14 . 1 1 the resulting biorthogonal functions are identical, 
but we will obtain different measures. If they are related using the flip of Proposition 14 . II we again have the 
same functions, but different measures, some for \q\ < 1 and others for |g| > 1. In Appendix \X\ we give a 
complete list of all different biorthogonal systems we obtain in this way and describe how they are related 
using these two kinds of symmetries. 



6. Measures 

In this section we briefly list the measures associated to vectors within the polytope P. The results in here 
are slight generalizations of the results in [2 , where we considered the limits of, essentially, (1,1). Plugging in 
functions does not significantly alter the arguments which shows that these limits are valid, as the functions 
we consider converge uniformly on compact sets outside some isolated points, by definition. 0] is devoted 
to extending this section to the multivariate case, and includes detailed proofs of the results. 

One important idea is to remember that there exist ruf and such that 

t, s r{uoz^^) r(uiz±i) 

I [Unq z^^) 1 (uiQ ™<ir±J-i 

are holomorphic. The bilinear form can then be expressed as 



(Lg) 



to-ti,t2,t3,ti,ta:q;p 



'^Uo<r<s<5^itrts;P,q) 

, lt=o nUz^^; p, q)T{uoq-"'fz^\uiq-^^z^^) dz 



f{z)g{z) 



c 
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r(z='=2;p, (7) 27riz' 



that is, the same measure with shghtly different parameters, multiphed by a holomorphic uniformly converg- 
ing (as p 0) function. For the integral limits below one will only see this idea reflected in the choice of 
contour, but for the series limits we need to actually consider these new functions. It is useful to introduce 
the notation 

ir:=tr, (0 < r < 3), := iio := q'^^f uq, := ui := q~"^^ui, 

for these new parameters. Moreover we introduce rrii for 0<i<5asmi = 0if0<i<3 and = ruf and 
m5 — nig. Notice that the lc{f) and lc{f) are related via 

lc{f){z) = lc{J){z){u^/z-q)]zy'''^UoZ-q)]zr'''' f f (" — )"^<z("')) ^^"^'^^ ■ 

' \ UQ J \ UqZ J 

Since / is holomorphic, we know that lc{f) is holomorphic, so this equation gives the possible locations of 
poles of lc{f). 

The first proposition is the analogue of ^ Proposition 4.1], and describes Nasrallah- Rahman type bilinear 
forms. 

Proposition 6.1. Choose generic parameters satisfying Ylr^r = q- Let a G R^, 'Y^^^o'^r = 1 cind ol^ > 
/or < r < 5. Let f G i(i4p"'') and g G A{t5p°"'). 
We now have the limit 



9) no<r<.<5:a.+a.=oM'*; 9) 

\JV',''rF JiyV'^'-rF ;/t,p=r— — — , +-U-1\ \ 

^Uo<r<s<5:a^+a^ = li<ltr ts )■. q) 



lim (/(^. trp"^),giz; trp"-))t^ 



X / Hf){z)Hg){z)- 



where lc{f) = lc{f{z;p'^^tr)) and likewise for lc{g). Here the integration contour C ~ C^^ is such that it 
includes the points qHr, (for < r < 5 with ar — Q and j > 0) and excludes their reciprocals. The contour 
can be taken to be the unit circle if \tr \ < 1 for all r with ar ~ 0. 

Notice that, while it is obvious that we have to consider the leading coefficient of / and g as the t^ change 
with p according to p°'''tr, we leave the z's invariant with p, i.e., we have C = 0. Note that we therefore 
exactly have that |C + 5I = min(i, ^ + a,-, ^ + a,- + ag + at), i.e. we evaluate the functions at the right C, 
given a (according to Theorem l5.2l) . In the following propositions we will have different values of C, but one 
can easily check that it is always the right value. 

By breaking the z — > 1/z symmetry in the integral before taking limits we can obtain limits for different 
values of a. This gives first of all (compare with [2 Proposition 4.2]). 

Proposition 6.2. Let tr £ C be generic, satisfying Yir^r = q- Let < a,b,c < 5 and let a G M^, and 

— i < C < satisfy X]r=o "^r = 1, eta + at + etc = C and ( < ar < — C /o*^ r = a,b,c and —( < < 1 + C 
for r ^ a,b, c. Let f G A{t4p°"^) and g G ^(^5^"^). Then we have 

p— !-0 

{a.b,c} {a.b,c} ,s^{a,b.c} {trts\q) 

Qr + Os— — 1 Qi- + as — 

Y{T.se{a,b,c}{qtr^t^^;q)Y{ '■^s {qtr^ts^;q) 

l\re{a.b,c}{q/Uz;q)l\r({a.b,c}{qz/tr;q) f (z^ q) Ilre{a,b,c}iqZ / tr] q)'^ lc = -i/2 

X / lc{f)izMg)iz)^^^ 



I\r£{a,b,c}{tr/z;q) JJr ^{a,b,c} 

(trz;q) 
dz 

X 0{qz/tatbtc;q): 



{qz^;q)l\re{a.b,c}itrz;q) 

\ a,=-l/2 



2TTiz 

where lc{f) = lc{f{p'^z)) and similarly for g. We also have the usual conditions on the integration contour 
(i.e. if {trZ;q) is in the denominator of the measure, we include poles of the form q'^/tr, k > 0, and if 
{tr/z; q) is in the denominator we exclude poles of the form tj.q"*' , k > 0). 
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The next case we consider is an integral as in [7, Proposition 4.4]. 

Proposition 6.3. Let tr E C be generic such that Yir^r — q- Let < a, 6 < 5 and let a G R^, and 

— i < C < satisfy X]r=o '^r — ^, eta = cub = C — C < c^r < 1 + C f^T r ^ a,b. Let f G A{t4p"'^) and 
g G i(i5P"0 

Then we have the limit (where C satisfies the usual conditions on a contour) 

p— »0 

_ (g; q)itatb; qy^ = -'/' Ur:a^ = -citrta,trtb.. g) 

Ylo<r<s<5iqtr^tJ^ ; q) 

a^^+as — 1 

X flc(f)(z)lc( )(Z) n.:.,.^l+c(g^A''-g) f 1-^^ y^^-^^^ ^(^^.df^;g) Cfz 

JC (ia/z,W^;g)nr:a. = -c(*'-^5«) V (*a2, tfe^; g) / e{taW,tbW;q) 2niz 

where w is some new parameter (which does not affect the value of the bilinear form when applied to our 
functions). Here lc[f) = lc{f{zp'^)) and likewise for lc[g). 

There also exists a double series expression for this bilinear form. The double series expression is the 
extensions to a measure of some of the cases in [21 Proposition 4.3]. It can be obtained from the integral 
expression by taking residues while shrinking the contour to the origin. In the case {a, &} n {4, 5} 7^ we 
get extra complications, as we would have to consider residues of lc{f) and lc{g). Thus our proposition 
statement will exclude those cases, though after the proposition we will explain how to get a more general 
expression. 

Proposition 6.4. Let tr E C be generic such that Ylr^r = g- Let < a, & < 3 and let a G K^, and 

— i < C < satisfy Y^^^o'^r — ^, eta = ctb = C ^.''^d, —( < ctr < 1 + C for r a,b (0 < r < 5). Let 
f G A{t4p"^) and g G i(i5p"^) 

Then we have the limit 

1 V^-i Ur:c.. = -citrtb;q)Ur:a.^l+A'-t'l) 



iqtl;q)J (|^;g)n.,.a.+„.=i(7^;g) 

.Y.lcm.g^MgXtaq'^) ( ^^^-'^I'^V^'^^'^^Y ^^'' . .."rrr"'^'^'"^'!.. 1' 



k>0 



(g^^;g)y (!r;g)a,.:„.+.,=i(t^;g) 



^i.r,-.n..M,.j.M...k, f {qtlqhk{tltbta;q)k\'^-i Ua^^-c(trtb;q)k 



X ^ lc{f){tbq^)lc{g){tbq^) --^'''-7.;;; Tn ^ 

where lc{f) = lc{f{zp'^)) and likewise for lc{g). 

To extend the result to the case a = 4 we should consider the series based around taq'' as a bilateral series 
which happens to have a term {q; q)k in the denominator which makes it zero for fc < 0. However, if a = 4 
we would be evaluating lc{f) at one of its poles (which are located at taq'' for —rrif < k < 0). As long as 
lc{f) indeed has a pole, it would cancel the zero from the {q; g)fc-term, and we would get an extra term. In 
particular we get at most mf extra terms. To make this more formal one could rewrite the summand in 
terms of lc{f) instead of lc{f), so that we can evaluate lc{f) in the relevant points. It should be noted that 
in this case the support of the measure becomes uoq^, not just itog^-", plus the support coming from the 
second series. 

The final limit case for the continuous measure corresponds to single sum expressions. Many of the issues 
surrounding the double sum expression arise arise here as well. Compare this proposition to the single series 
cases of [21 Proposition 4.3]. 
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Proposition 6.5. Let tr E C be generic such that Yl,. tr = Q- Let X]r=o Q^r = 1- Let < a < 3 be such that 
^\ < «a < and 1 + Ola > otr > Cia foT r^a, 0<r<5 and such that 1 > + > for r,sj^a. 
Moreover assume 2aa = 'l2r^a-ar-+aa<oi'^r + eta)- Let f G A{t4p°'*) and likewise g € A{t5p°'^). 



We have 



limp— ^/(^. i^pa.), ^(^. i^pa.)) 
p— >0 



■\~tlg^^ 



1 1 



X 

fe>0 



^Zc(/)(t„g'=)k(<?)(i.g'=) [—^{tl-^q) 



a/tr,q)k ^ ' 




where N =\\r \ r ^ a,ar < —aa}\, and lc{f) = lc{f {zp°''^)) , and likewise for lc{g) and the valuations. 

Finally let us consider finite measures. In this limit we want to insist that toil = for all p, so 

ao + ai = 0. The proof in this case is just interchanging sum and limit. Note that the resulting expression 

looks very much like the sum cases above. 

Proposition 6.6. Suppose tr E C are generic, such that toti = q~^ and t2t2,tAt5 = q^~^^- Suppose 
— 5 < ckq = —0^1 < 0, and 1 + ao > oir > cxo for r > 1, and that + Q!s < 1 (for 2 < r < s < 5) and 
Q!2 + as + a4 + as = 1 and J22<r<5:ar<-ao ao + ar = 2ao- Let f e A^tip""*) and likewise g e A^t^p"^). 
Then (/, g)t,u € M{i, u, q, t) and 

N 

lim p— ))— )) (/, = V lc{f) {toq^Mg) {toq')wk,a {U), 
«— >o ^ — ' 

fe=0 

with lc{f) = lc{f{zp°'°)) and lc{g) = lc{g{zp°'°)) . Here the weights Wk,a are given by 
• //ao = 

_ ^-^og^M9~^:io;g)fc / 1 ^^-2^ 1 TT {totr;q)k{titr;q)N , ^fcg^S) 

l-*o {q,qto/ti;q)k\titlqJ 2<r<57c«,=o il*o/tr;q)k ° 

{totr;q)k{tltr;q)N , , , x_fc. , , .-N -C'^-C^) TT 1 



X 

2<r<5:c«r=l 



(qto/tr; q)k 2<.<.<5:t+a3 = l ^"^'^^^^'^ 

7/-l/2<ao<0 

TT MVffl., , . TT I ,,^k (-) 



{q;q)ktW^'^) 2<.<il=„o 



2<r<b:(XQ<otr<—oto 



X n M^;^)^ n n 



2<.<5 a.=-ao 2<r<5 i.=l+ao (l*o/*r;q)k 2<,<,<il+„, = l 5)iV 

• If ao = -1/2 

^^^^ l-tWHq-'^,tU). 1 N±V {J±\\.i^).0 

1-^0 {q,qto/ti;q)k {ti/to;q)N \ti J V 

(io^r; q)k{titr; q)N{-qtQ/tr)''q^^^ 



X 



, [qto/tr] q)k 



1 

X 



2<r<5:a.=l/2 (^^oA- 2<.<«<5 ^^A'^*- 
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7. The g-AsKEY scheme 

In this section we want to highlight the part of our degeneration scheme of biorthogonal famihes which 
corresponds to the g-Askey scheme. The g-Askey scheme is a unified way of presenting the classical q- 
hypergeometric orthogonal families of polynomials, and put them in a picture which clarifies the possible 
limit transitions. Orthogonal polynomials are just the special case of biorthogonal rational functions, where 
the poles of the rational functions have disappeared to infinity and the a priori two different families on both 
sides of the bilinear form have reduced to identical functions. Therefore it comes as no surprise that the 
different families in the g-Askey scheme return as special cases in our degeneration scheme of biorthogonal 
families of rational g-hypergeometric functions. 

We would like to emphasize however that we have a slightly different philosophy on what we include in 
our scheme and what we consider to be different. In particular we do not care about positivity of either 
the squared norms, or the measure, whereas the g-Askey scheme does. It only makes sense to worry about 
this in the context of orthogonality (as opposed to biorthogonality) , but in the case of orthogonality it is a 
very important property. Because of this difference we include families which have no positive measure at 
all, while the g-Askey scheme would consider two families of orthogonal polynomials to be different if they 
are the same functions, but in different ranges of the parameters, for which the conditions of positivity are 
different. This always happens if different measures work for \q\ < 1 and |g| > 1 (e.g. compare Al-Salam 
Carlitz I and II), which in our scheme means that a degeneration is not equal to the degeneration of the 
flipped version. However it also happens for just simply different parameters while keeping |g| < 1, for 
example, consider the relation {q;q)nLi"\—q^^]q) — C„(g~^; — g~"; g) between g-Laguerre and g-Charlier 
polynomials. 

A special case of this is when we can specialize the product of two parameters to g~^ which makes 
the measure finitely supported, for example compare Askey- Wilson polynomials and g-Racah polynomials. 
These cases will be the same point in our degeneration scheme. However it is relatively easy to determine if 
we can specialize parameters in this way. This is only possible if we have two visible parameters ar which 
sum up to (because their product has to be g^^ for all values of p). 

The degeneration scheme in Figure [1] shows the different limits we obtain of families of orthogonal poly- 
nomials. We modded out by the permutation symmetries in the parameters on the elliptic level, and the 
translations along p-shifts. However we did not mod out by the flip (sending g — > 1/g), unlike in the full 
scheme in Figure [21 We also changed the names from that scheme. The iterated limit property shows that 
all families are limits of the top level, the Askey- Wilson polynomials. For each degeneration there is a unique 
associated face which contains {a; J'X) — ^ ■— (0000; ^i; 0) at its boundary. This vector corresponds to the 
Askey- Wilson polynomials themselves. In the midpoint of the face denoted by [a/b] we have 70 = 7i = a/b. 
Moreover the degeneration [a/b] is a (5 — 4)/2 level degeneration down from the Askey- Wilson level. The 
level of degeneration and the value of 70 = 71 almost completely determine which degeneration we are in, 
except that there are two limits corresponding to [5/10], so we called one of them [5/10]'. 

We can also put the different degenerations in a list, giving their identification with elements of the usual 
g-Askey scheme. The table lists the names we gave the degeneration (both in our scheme of orthogonal 
polynomials, and the complete scheme of biorthogonal functions), the midpoint of the face giving this limit 
which contains v, and the classical polynomials it corresponds to (both for measures with infinite support 
and, if it exists, with finite support). These families can all be found in [6], except for the ones on the right 
most line (below the top level), and the ones for which no positive measure is possible (denoted by NP in 
the scheme below). The case [6/8], i.e. Al-Salam Chihara polynomials with jgj > 1, was studied by Askey 
and Ismail [1]. 
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[2/10] 



[2/12] 




[8/10] 



[10/12] 



Figure 1. The g-Askey scheme as part of the degeneration scheme for biorthogonal q- 
hypergeometric rational functions. 
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8. Pastro Polynomials 

We would like to draw attention to the special points in the degeneration scheme for which the lim- 
its are families of biorthogonal polynomials (i.e. polynomials outside the g-Askey scheme where we have 
orthogonality). These polynomials are the Pastro polynomials 8 , and its degenerations. 

The top level of these polynomials is llllpp, associated to the vector a = {—j,0, j, ^;0, ^) (which equals 

it's own flip). For the first function i?„(zp^3 ; toP^^ , ii, ^2^^ , isP^; "o, "iP^) the valuation turns out to be 
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zero, and we get the explicit expression 

(9) P„(z) := lim Rn{zp^i;toP^Kti,t2pKt3P^;uo,uipi) 

The 3(^2 expression is obtained by interchanging hmit and sum in the definition ^ of the biorthogonal 
function The 20i expression then fohows using [5] (III. 7)]. 

Notice that the parameters of the function only appear in certain combinations. In particular if we define 

Pn{w; A, B) ^ —— — 201 9-n/n 5 9 



we can write 

1 /2 

q ' q q 

Pn{z;to ■■ ti,t2,h]UQ,Ui) =Pn{— ; T , T )• 

t2huiz tiuo tsui 

For the right hand family (the functions with uq and ui interchanged) we notice that by the symmetries 
from Lemma |3.2[ and the permutation symmetry in the t^'s we have 

Rn{zp~i;top~i , ti,t2P^ , hp^ ; uip^ , uq) = Rn{zp^ ; t^p^ , tip^ , t2P~^ , h] ui,uqp^) 

= Rn{z^^p^^ ; toP^ , tiP^ Mp~^ ui,uop^) 

e{p-^t2h,P^tit2,pi^,P~TjZ:-)n~ 1 _i ^ _i _ 1 ^ 1 1, 

= —-^ 1 Zr~i 3ot Rn{Z P ■^]t2P 'i,t3,top-i,tip^;Ui,UQP^). 

9{p4toti,p^tot3,p ^J^,P*^)n 

As a corollary we see that the limit on the right hand side is up to a constant equal to the limit on the left 
hand side with different parameters. Indeed we have 

1 1 111 / 1 \ " 1 

Qn := \imRn{zp^^;top^^,ti,t2P^,t3P^;uip^,uo) = — I P„(-; t2, ^3, ^o, ^i; mi, wq). 

p^Q \Uotot1t2 J z 

In particular, setting 

qniw;A,B) :=p„(-;B,A), 

w 

we have 

/ (7 \ " q^^^ q q 

Qn{z;to : ii,t2,t3;uo,-Ui) = ; t , t )• 

\huij t2huiz tiuo hui 

Note that this is the same parameter correspondence as between Pn and Pn- Direct comparison now shows 
that the univariate functions Pn and g„ are equal to the polynomials [8j (3.1)], up to a constant and under 
the parameter correspondence q°' — A, q^ = B and t — w. 
The limit of the biorthogonality relation gives us 

(P„(-; tr\Ur\ q), Qm{-; U; Ur, q)) = (5„.m(ti'«o)""7;%^^^- 

{tot2', q)n 

where 

,f V {q,tot2;q) [ ft \ I J{qz/tQtiUQ]q) dz 

{q/tiUo,q/t3Ui;q) Jc [to/ z,t2z; q) 2mz 

or in terms of p„ and q.n this becomes 

{p^{.-A,B),q^{--A,B))=5n,m(—\ 7^%, 

where 

, (g.^jg) e{q'/^w;q) dw 
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We would like to consider a special case. Indeed we want to specialize t^ui 1, or equivalently B ^ q. 
The simplest expression to do this in is in the original 3^2 of (l9|), which then becomes a 2^1 which can be 
smnmed by the g-Vandermonde j5j (11.6)], to give 

that is, the polynomials p„ become just the monomials w". Observe that the specialization A^B ^ q turns 
the limiting measure into the trivial , and indeed we know that the orthogonal functions with respect to 
this measure are and w"". The multivariate analogues of this special case of the Pastro polynomials are 
the Macdonald polynomials 7] (indeed the univariate part of the measure associated to the orthogonality of 
the Macdonald polynomials is the trivial dw/2'Kiw). 

Appendix A. Description of the limiting systems 

In this section we write down explicitly all the limiting families of biorthogonal systems, and try and 
give a universal picture of the different cases. After the theory of Section [5] this has reduced to a concrete 
combinatorial problem. 

In Section [5] we have seen that the limits of pairs of biorthogonal functions, together with measure, are 
classified by the different kinds of faces of the splitting of polytope P^°^ in Pj, Pn and Pm (as in ,2 ) 
(with C given by its prescribed value), with the exception of the interior of P//, modulo the symmetry action 
54 X 5*2. Moreover we also still want to mod out by the shifting symmetries (while we are looking in a 
fundamental domain for this shift-action, boundary points can of course still be mapped to other boundary 
points). On the other hand we consider two solutions associated to each other by the flip as different (as the 
measures only make sense if we impose \q\ < 1, while the flip maps g — > 1/q, thus making it impossible for 
this condition to be satisfied for both cases). 

The different kind of limits are named according to the relation of the a to C. If ^ = or ^ = 1/2 we have 
two integers in front, the first one counting how many (0 < r < 3) are equal to C (or some integer shift 
of (), the second counting how many are unequal to ( (and its integer shifts). If C 7^ and ( ^ 1/2, we have 
four integers. The first two counting how many are equal to ( (and its integer shifts), respectively — C (and 
its integer shifts). The second two counting how many are in the interval (— C,C) (and shifts), respectively 
(C, 1 — C) (and shifts). These two pairs of integers are sorted large to small. The letters at the end describe 
the relation of 70 = a4 and 71 = as to C- If (modulo 1) a4 = = ±C then we write v2, if 04 = —as = ±C 
we write vv, if q;4 = ±C and as ^ ±C or vice versa, we write vp, if a4, as G (— C, C) or a^, as G (Ci 1 ~ C) 
(and shifts), we write os, and if a4 G (— C, C) and as G (C, 1 — C) or vice versa we write pp. 

These names give a rough description of the associated biorthogonal functions. Indeed if a^ = ±C 
(0 < r < 3) then we will be able to "see" this parameter in the limit, otherwise the limiting function will 
not depend on it. If two a^'s are equal the limit will moreover be symmetric under their interchange. If 
a4 = ±C or as = ±C then the associated biorthogonal functions will be rational functions, otherwise they 
will be polynomials. Again if they are equal or in the same interval, as in the cases v2 and as, the two 
families will be related by interchanging uq and ui. The pp limits are Pastro Polynomials and their limits, 
while the as limits are limits in the Askey-Scheme of g-hypergeometric orthogonal polynomials. 

The simplest way to arrive at this classification is to start with the different faces mod Sq of Pi , P// and 
Pjii from [2], and tally their 5*4 x ^2 orbits. Subsequently we identify those faces which differ by a shift. 
The flip symmetry is treated separately, as it involves inverting q and thus alters the feel of the functions. 
All faces are simplicial, so we will just give the different vertices of the faces, together with the midpoints. 
Finally we briefly indicate the type of measure, whether it is a Nasrallah-Rahman type beta integral (NR) 
(Proposition l6. 1| ) ■ a symmetry broken integral (SB) (Proposition l6.2| ). a single series measure E (Proposition 
16. Sp . or a double series measure (Y?) (Proposition 16.31 these double series can also have an expression as a 
symmetry broken integral). 

We will use the following notations for the vertices: dj {0 < j < 3) is the vertex with all O's except aj = 1 
(so IC+ 5I — 5)- {j = 0, 1) is the vertex with all O's except = 1 (with again |C+ ^1 = 5)- fij is the vector 
with all ^'s, except a^ = a^ = — ^ (and |^ + ^| = 0). gij is the vector with all ^'s, except a^ = 7j = — ^. 
Finally hoi = (i, i, i, i; -i, -i; 0). 

A.l. Top level. 
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A. 1.1. 40v2. At the top level we consider vertices of Pj, Pjj, and Pju and get up to ^4 x ^2 symmetry the 
vectors ^3, ei, /oi, 500 and hoi- These vectors are all mapped to each other by purely shifts, so in fact we 
only have one case. However the different vertices do correspond to different measures with respect to which 
the limiting rational functions are biorthogonal. 



A. 2. Level 2. Let us tabulate the different edges up to the S4 x S2 symmetry. Then we group those which 
are related to each other by shifts in the W{Ee) lattice. The midpoints of the faces are given as the vector 

(ao, ai, a2, «3; 70i 71; ^0- (Wc negated ( as in the polytope C is always positive, but by the z — > 1/z 
symmetry of the elliptic hypergeometric biorthogonal functions this negative value of ( is just a consequence 
of some choices we made). 



Name 


Vertices 


Midpoint 




31vp 


ei,d3 
foi , 500 
5007 ^01 


(0,0,0.i;0,i;0) 
(-1 I 1-0 1-1) 

\ 2' '2'2' '2'2'' 

V"' 2 ' 2 ' 2 ' 2 ' ' 2 / 


NR 
E 
S 


2200i;p 


ei,/oi 


/ 1 111.13.1^) 

V 41 41414:4:4:4/ 
/ 1113. 11. l^ 

I 4'4'4'4' 4'4'4/ 


S2 

S2 


3100t;2 


61,500 

d3,foi 

da, hoi 


/ 1111. 13. 1\ 

I 4'4'4'4' 4'4'4/ 
( 1 113.11.1\ 
I 4' 4' 4' 4' 4' 4' i) 
^1113. 1 1.1\ 

\4'4'4'4' 4' 4'4'' 


S2 
S2 
S2 


22v2 


d2,d3 
foi , fo2 


(0,0,i,i;0,0;0) 

(-_! i-i) 

V 2' ' '2'2'2'2-' 


NR 
E 


AOas 


eo,ei 
500,501 


(0,0,0,0;i,^;0) 
f-i 1 1 i-0 0- 

V 2' 2' 2' 2' ' ' 2/ 


NR 
E 



The case 40as is the top level case in the 5-Askcy-schcme and corresponds to Koornwinder polynomials, or 
univariately, Askey- Wilson polynomials. We obtain two different measures, the one associated to the edge 
between eo and ei is the usual Askey- Wilson measure. 



A. 3. Level 3. Now wc consider 2d faces, that is, triangles. In this case we have some faces which map to 
different faces if wc perform the g — )■ 1/g flip. Considering they are essentially the same functions, but also 
realizing that we consider them on different domains (either \q\ < 1 or \q\ > 1) we group them together but 
are aware of the differences. 
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Name 


vertices 


Ordinary 
midpoint 




vertices 


Flipped 

midpoint 




22vp 


^2,^3, ei 
./oi- .'7(10- .'710 


(0,0.i,i;0,i;0) 

( 1 1 ' 1 ' 1 . ' 1 1 . 1 \ 


NR 
SB 


.foi, .fo2, .900 
,'/()()■ .910- /'nl 


/ 1111. 11. In 

V 2' 6' 6' 2' 6' 2' 2/ 

(1111. 1 1 . 1 \ 


s 


V (i ■ () ■ 2 ■ 2 ■ (>■ 2' 2) 


W) ■ () ■ 2 ■ 2 ■ 2 ■ (i ■ 2 ' 


2110vp 


t^3, ei, /oi 

C^Sj .900; ^'01 

ei: 500j /oi 


( 1 1 1 1 . 1 1 . 1 \ 
^ 6' 6'6'2'6'2'6'' 

v"' 3' 3' 3' 3' ' 3' 

V 3 ' ' 3 ' 3 ' ' 3 ' 3-* 


Y? 
S 
S 


£^3,61,, 900 
•^3, /oi,5oo 


/ 1111. 1 1 . 1 N 
6 ■ 6 ■ 6 ' 2 ' 6 ' 2 ' 6 

f-i i 2-0 i- i) 

V 3' ' 3' 3' ' 3' 3/ 


E 


lUOvv 


(^2,^3,500 
'^3 7 500,fflO 

ei, /oi, /o2 


/ 1111. 1 1 . 1 \ 
^ 6'6'2'2' &^ & ■ g' 
/Q 1 2. 11. In 
V"' 3' 3' ?' 8' 

C-i i -• i") 

V 3 ' 3 ' 3 ' 3 ' 3/ 


Y? 
S 
S 








2020v2 


d2,d3, foi 

ei,ffoo,5io 


/ 1 1 1 1 . 1 1 . 1 \ 

^ 6' 6'2'2'6'6'6'' 
/n Q 1 1. 1 2. In 
v"' 3' 3' 3' 3' 3/ 




£^2, C?3, 'iOl 

c'a, /oi, /02 


/llll. 1 l.lN 

V6'6'2'2' 6' 6^ 0' 

(—l {) 1- 1 1- 1\ 

\ 3 ' ' ' 3 ' 3 ' 3' 3'' 


S 


13v2 


di,d2,rf3 

/oi, /02, /l2 


(0,i,i,i;0,0;0) 

/ 1 1 

V 6' 6' 6' 2' 2' 2' 2/ 


NR 
SB 


/oi, /02, /o3 

500,510,520 


( 1111. 11. In 
V 2'6'6'6'2'2'2/ 

/I 1 1 1, 1 l.lN 

^6' 6' 6' 2' 2' 2' 2/ 


S 

s 


31as 


c?3,ei,eo 
.900, .901, hai 


(0,0.0,i;i,i;0) 

/ 1111. 1 1 . 1 N 

V 6'2'2'2' 6' 6 '2-' 


NR 
SB 


/oi, 500, 501 


/ 1 1 1 1. 1 1. In 

^ 2' 6' 2' 2' 6' 6' 2) 


s 


2200as 


ds, 900,901 
Co, ei, /oi 


V 3 ' 3 ' 3 ' 3 ' ' ' 3-* 

1 1 lll.ll.lN 

N 6' 6'6'6'2'2'6'' 


S 

S2 









The points where we do not include flipped versions, in this table or the coming ones, are all their own 
flips. In particular for those points the same pair of families of functions is biorthogonal with respect to 
explicit measures both when jgl > 1 and when \q\ < 1. 

A. 4. Level 4. The 3-d faces are all tetrahedra. 



Name 


vertices 


Ordinary 

midpoint 




vertices 


Flipped 

midpoint 


M 


13vp 


di,d2,d3,ei 




NR 


foi, f02, /03, 500 

900,910,920, hoi 


( 1 1 1 1 . 1 1 . 1 N 

I 2'4'4'4'4'2'2'' 
(-1111. 11. In 
V4'4>4>2' 2'4'2-' 


Y 


2020?;p 


d2, ds, ei, /oi 
ei,/oi,5oo,5io 


/ 1 1 3 3. 1 3. In 

^ 8' 8'8'8'8'8'8'' 
( 1 13 3. 1 5 . 3 N 
V 8' 8'8'8' 8'8'8'' 


T? 
SB 


da, foi, f 02,900 


C 3 1 1 5.1 3. 3n 
V 8'8'8'8'8'8'8^ 


Y 


1021t;}j 


c^2,d3,/oi,5oo 

c?3,ei,5oo,5io 


V 4i'-')2'2' '4'4'' 

CO 1 1- 11- 1) 


s 


d2, da, goo, hoi 
d3,ei, foi, fo2 


CO i i i--i 0-i) 
C-i i- i i- i) 

V 4 1 2 ' 4 ' 2 ' 47 


E 

E 


lOBOvv 


di,d2, ds, 900 
ei, /oi, .fo2, /12 


( 1 3 3 3. 1 1 . In 
V 8'8'8'8' 8'8'8^ 
/ 1 1 1 3 . 3 5 . 3 N 

V 8' 8' 8'8'8'8'8^ 


Y? 
SB 


da, 900,910,920 
ei, /oi, /02, /o3 


/I 1 1 5. 3 3. 3n 
^8'8'8'8' 8' 8' 8' 
( 3 1 1 1.3 5. 3N 
V 8'8'8'8'8'8'8'' 


E 
E 


llllpp 


ds, ei, /oi, 500 


C-i 1 i-0 I- 1) 

V 4''4'2''2'4/ 


S 








n20vp 


^2,^3,61, 500 

c^3,/oi,5oo,5io 


( 113 3. 1 3 . 1 N 
^ 8'8'8'8' 8'8'8-' 
/ 1 13 5. 1 3. 3n 
V 8' 8'8'8' S'S'S/ 


SB 


'^3,5oo,5io,/ioi 
ei, foi, f 02, 900 


(1 1 3 5 . 3 1 . 3 N 
N8'8'8'8' 8'8'8'' 
( 3 1 1 3. 1 5. 3n 
V 8' 8' 8' 8' 8' 8' 8/ 


E 
E 


0022t;w 


d2, ^3,500,510 


CO 1 1- 11- 1) 

u,2,2, 4'4'4/ 


s 








1030t;2 


dl, ^2, C?3, ^01 

da, foi, f 02, fi2 


(-1 3 3 3.1 1. In 

^8'8'8'8' 8' 8'8^ 
/ 1 1 1 5. 3 3. 3 N 
^ 8' 8' 8'8'8'8'8'' 


SB 


ei, 900, 910,920 


(-1113. 3 5. 3n 
V8' 8' 8' 8' 8' 8' 8/ 


E 


04v2 


dQ,di,d2, ds 


Ci i i i-0 0-0) 

V 4 ' 4 ' 4 ' 4 ' "i 


NR 


500,510,520,530 


/-I 111. 11. IN 
1,4, 4, 4, 4, 2'2'2'' 


E 


22as 


d2,d3,eQ,ei 


CO i i- i i-0) 


NR 


/oi, /o2,5oo,5oi 


( 1111. 11. In 

V 2'4'4'2'4'4'2' 


E 


2110as 


da, eo, ei, /oi 
c^3,5oo,5oi,/ioi 


/ 1 1 1 3. 3 3. In 

V 8' 8'8'8'8'8'8'' 
1 3 3 5. 1 1 . 3 N 

^ 8'8'8'8' 8' 8' 8/ 


SB 


t^3,/oi,5oo,5oi 


I 3 1 3 5.1 1.3N 

V 8' 8' 8' 8' 8' 8' 8/ 


E 


1120as 


d2,d?,,9oo,9oi 
eo, ei, /oi, /02 


C-i i i i-0 0-i) 

I 4, 4, 2 ' 2 ' 47 
C_i i-i i-i) 

I, 4)"-')")4)2'2'47 


S 

s 
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A. 5. Level 5. The 4-d faces are all simplices as well. To save some horizontal space we use an abbreviated 
notation for the midpoints. For example (| ; 0, -g; 0) denotes the point (-g, -g, -g, -g; 0, 1; 0). 



Name 


vertices 


Urdmary 

midpoint 


M 


] 

vertices 


^ lipped 

midpoint 




OAvp 


dn d^ d'? d-i e^ 


V 5 , 5 , 


NR 


(7nn (7i n (79n 0'?n /7,ni 


(-34. 1 3 . i\ 

V 10 ' 2' 10' 2' 


T, 


OOSlwp 


di d'? d'i Onn hm 
<^3i ei, /oi, /o2, /l2 


(0 23._1 0;1) 

V"' 5 ' 5 ' ' 5' 
/ 1 3 1 . 3 1. 3 \ 
V 10 ' 2' 10' 2' 10' 


SB 
SB 


d'i ei Onn Oi n 09n 


/ 1 3 1 . 3 1. 3 \ 
V 10 '2' 10' 2' 10' 


E 


1021pp 


do d'i e^ fni Onn 

^o; '-'I7 i/^J^J 

rf3,ei,/oi,5oo,ffio 


V 5 ' ' 5 ' "' 5 ' 5 ' 
( 1231. 1 1 . 3 ^) 
\ 10 ' 10 ' 2 ' 10 ' 2 ' 10 ' 


SB 
SB 


fii r/^ /ni /no Onn 


( 3 1 2 1 , 1 1 . 3 \ 

\ 10 ' 10 ' 2 ' 10 ' 2 ' 10 ' 




0022t;p 


(^2, 0^3,61, 300) 510 

d2,d3, /oi, ffoO) 510 


f02 22. _1 2.U 

' 5 ' 5 ' 5 ' 5 ' 
/ 1 2 l2^ 1 3 . 3 \ 
V 10 ' 2 ' 10' 10' 10' 


SB 
SB 


C^2,C^3,500,5l0, /iQl 


(- 1 2 l2, 3 1 . 3 ^ 

\ 10 '2 ' 10 ' 10 ' 10 ' 


S 


0040i;2 


do,di,d2,d3, hoi 


^34. 1 2, 1 . 
V 10 ' 10 ' 10' 




ei, 500,510,520, 530 


/l-i. 2 3. 2n 

V5 ' 5' 5' 5' 


S 


W30vp 


di,d2, dz, ei, 500 


( 1 33. 1 3 . 1 \ 

^ 10' 10 ' 10' 10' 10' 


S2 


'^3,500,510,520, /loi 
ei, /oi, /02, /o3, 5oo 


(1^ 3. 2 1. 2\ 

'5', 5' 5' S-' 
/ 2 l3. 1 3. 2\ 
y 5' 5 ' 5' 5' 5' 


s 


13as 


^1,^2,^3,60,61 




NR 


/oi, /02, /o3, 500,501 


f 1 33.321- 
V 2 ' 10 ' 10 ' 2 ' 




2020as 


d2,d3, eo, ei, /oi 


( 1 2 3 2 3 2 2 N 

^ 10 ' 10 ' 10 ' 10' 


S2 


C^3, /oi, /02,500,501 


/ 2123. l2. 2\ 

y 5 ' 5 ' 5 ' 5 ' 5 ' 


E 


1021as 


o?3, eo, ei, /oi, /02 

^2,^3, 900,501,^101 


f-l 0^ 2.22.1^ 

V 5 ' ' 5 ' 5 ' 5' 
( 1 3 l2. 1 2. 3 X 
V 10' 10' 2 ' 10 ' 10' 


SB 
SB 


^2,^3,701, 500,. 901 


( 3 1 l2. 1 2. 3 s 

V 10' 10' 2 ' 10 ' 10' 


E 


lOaOas 


rfl,'^2,rf3,500,501 

Co, ei, /oi, /02, /12 


f-l 23.02. l^ 

V 5' 5 ' 5' 
/ 1 3 3 . i2, 3 N 

V 10 ' 10' 2 ' 10' 


SB 
SB 


Co, ei, /oi, /02, /o3 


/ 3 1 3. i2_ 3 s 

^ 10' 10 ' 2 ' 10' 


E 



A. 6. Level 6. The 5-d faces are not all simplices anymore, but the only non-simplices are the interiors of 
P// t, which are excluded anyway. 



Name 


vertices 


midpoint 




0022pp 


'^2,c?3,ei,/oi,5oo,5io 


/ 1 15 5. 1 5 . 1\ 

V 12' 12' 12' 12' 12' 12' 4-* 


SB 


OAas 


(io,rfi,rf2,d3,eo,ei 


(1 1 1 1.1 l.Q) 

'-e'e'e'e'e'e'' 


NR 


ooaias 


di,d2,d3, 500, 5oi, hai 
d3,eo, ei, /oi, /02, /12 


(-1 5 5 5. 1 1 • 1^ 

V 12' 12' 12' 12' 12' 12' 4' 
( 1 1 1 5 . 5 5 . 1\ 
12' 12' 12' 12' 12' 12' 4' 


SB 
SB 



In these cases the flips of the points are not in the polytope P*-"-* (thus in particular these points are not 
invariant under the flip). This implies that for the pairs of families of polynomials under consideration the 
measures from [2] only work when |g| < 1, not when \q\ > 1. In [4j we show that we can give measures in 
these cases also for |g| > 1, the measures then being related to bilateral series. 

A. 7. Degeneration scheme. Having tabulated all these different families of biorthogonal functions it be- 
comes possible to draw them all in a single graph, see Figure [2] 

The graph should be read as follows. In the graph we placed all families of biorthogonal functions except 
for the as ones (for those: see below). The height of a family denotes the dimension of the face which 
the family corresponds to, the top level (40u2) corresponds to a vertex, so dimension 0, while the bottom 
level (0022pp) corresponds to the interior of a Pm, which has dimension 5. This dimension is inversely 
proportional to the number of free parameters the limiting family still has. Indeed the top degeneration 
depends on 5 — = 5 free parameters (remember we have a balancing condition, which defines the sixth 
parameters in terms of the other five), while the bottom degeneration has no (5 — 5 = 0) parameters. 

The edges between different families denote limit transitions. By the iterated limit property we know that 
we can take further limits of our biorthogonal families. One family is a limit of another one if one of the 
associated faces of the first family contains a face associated to the second family. This allows us to easily 
obtain the possible limit transitions. 
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(22w2) 



(13i;2) 



(04w2) 




(l021pp) 1030UJ9 0022wp 



0040w2 




Figure 2 . The degeneration scheme of biorthogonal systems 

We indicated the Pastro-polynomial type hmits (i.e. the ones of form pp) with an ovaL The rectangles 
actuahy contain much more subtle information. Underneath (i.e. one level lower) each of the limits denoted 
with a rectangle there exists a family with the same node, except that the u parameters are of the as type. 
Moreover there exists a limit from the family in the box to the corresponding as- type limit, and this is the 
only limit from a non-as type family. The limits between different as-type families are given by the edges 
between the corresponding families with a rectangle around them. 
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